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In some recent articles we developed a new systematic approach to generate solvable rational 
extensions of primary translationally shape invariant potentials. In this generalized SUSY QM 
partnership, the DBT are built on the excited states Riccati-Schrodinger (RS) functions regularized 
via specific discrete symmetries of the considered potential. In the present paper, we prove that this 
scheme can be extended in a multistep formulation. Applying this scheme to the isotonic oscillator, 
we obtain new towers of regular rational extensions of this potential which are strictly isospectral to 
it. We give explicit expressions for their eigenstates which are associated to the recently discovered 
exceptional Laguerre polynomials and show explicitely that these extensions inherit of the shape 
invariance properties of the original potential. 
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I. INTRODUCTION 



Since the pionnering work of Gomez-Ullate et al jT| , the exceptional orthogonal polynomials (EOP) and their con- 
nexion with rational extensions of solvable quantum potentials have been an active research subject [2M20I] . The EOP 
Oh 1 appear to be the constitutive elements of the eigenstates of solvable rational extensions of the second category primary 
translationally shape-invariant potentials (TSIP) [2lT[24j In a series of recent papers [l8l - [20| we have developped a 
new approach which allows to generate infinite sets of such regular solvable extensions starting from every TSIP in a 
very direct and systematic way without taking recourse to any ansatz or ad hoc deforming functions. This approach 
is based on a generalization of the usual SUSY partnership built from excited states. The corresponding Darboux- 
Backlund Transformations (DBT), which are covariance transformations for the class of Riccati-Schrodinger (RS) 
\ equations [21( , are based on regularized RS functions corresponding to unphysical (negative energy) eigenfunctions of 
■ the initial hamiltonian. They are obtained by using discrete symmetries acting on the parameters of the considered 
family of potentials. If the use of negative energy states has been already proposed in early years of SUSY MQ 
development (25|, such a systematic scheme has never been envisaged. In the case of the isotonic oscillator, this gives 
the three infinite sets LI, L2 and L3 of regular rationally solvable extensions of this potential, the LI and L2 series 
potentials being strictly isospectral to the isotonic one and inherits of its shape-invariance property. Very recently 
Gomez-Ullate, Kamran and Milson |26( introduce a new family of EOP built from the classical Laguerre polynomials 
via two step Darboux transformations. In the present article we show that the scheme developed in |18l - l2G | admits a 
very natural n-step extension which allows to recover and generalize the result of [26| . We obtain new infinite towers 
of rational extensions of the isotonic potential strictly isospectral to this last. Their eigenfunctions are associated to 
J> . the new EOP families, the one introduced by Gomez-Ullate et al corresponding to the n — 2 particular case. 

V , 

The paper is organized as follows. First we recall the basic elements of the method in its one step form and the 
results obtained for the isotonic system. The regularity of the RS functions is here obtained in a self consistent 
way, using disconjugacy properties of the Schrodinger equation [27i - [29| without recourse of the Kienast-Lawton-Hahn 
theorem |30l.l3lT|. Then we present on a formal level the general m— step scheme and in the specific case of the isotonic 
oscillator, we give a criterion of concrete application in terms of asymptotic behaviour near the origin. Next we 
consider in a detailed manner the particular m = 2 case, initially envisaged by Gomez-Ullate et al (26|. We then prove 
by induction a sufficient condition to be able to build a to— step chain of regular and strictly isospectral extended 
potentials. We finally show that the shape invariance properties of the isotonic potential are hereditary and reached 
by all these extended potentials. 

During the writing of this article, two papers on the same subject appeared. The first one [32| contains results 
parallels to those developed here. The authors adopt a point of view similar to the one initiated in [l8-20] and 
extended in the present paper. The so-called "virtual states deletion method" corresponding exactly to the use 
of DBT based on unphysical eigenfunctions associated to regularized RS functions. The second one [33| discusses 
quite extensively the 2 and 3 steps cases under a slightly different angle and proposes an interesting conjecture on 
possible "degeneracies" in the set of the new generated ELP. The content of the present article has been established 
independently of these two works. 
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II. DARBOUX-BACKLUND TRANSFORMATIONS (DBT) AND DISCRETE SYMMETRIES 



If ip\(x; a) is an eigenstate of H(a) = —d 2 /dx 2 + V(x; a), a £ R m , x £ I C R, associated to the eigenvalue E\(a) 
(E o (a) = 0) 

^(x; a) + (E x (a) - V(x; a)) ^(x; a) = 0, (1) 

then the Riccati-Schrodinger (RS) function w\(x;a) = —ip'^(x;a)/rp\(x;a) satisfies the corresponding Riccati- 
Schrodinger (RS) equation [21| 

- w' x (x; a) + wl(x; a) = V(x; a) - E\(a). (2) 

It is a well-known fact that the set of general Riccati equations is invariant under the group Q of smooth SL(2, R)- valued 
curves Map(R, SL(2, R)) [UllEl- The particular subclass of Riccati-Schrodinger equations is, as for it, preserved by 
a specific subset of Q. These transformations, called Darboux-Backlund Transformations (DBT), are built from any 
solution w 1/ (x; a) of the initial RS equation Eq© as [2l|, [34| HH 



/ s A(w v ) (y). . / \ i E x (a) - E v {a) 

w\(x;a) -» w\ '{x; a) = —w v (x; a) H r r, (3) 

w u (x; a) - w\(x; a) 

where E\ (a) > E v (a) . is then a solution of the RS equation: 

- w^'(x; a) + (w{ u) (x; a)) * = V™ (x; a) - E x (a) , (4) 
with the same energy E\(a) as in Eq([2]) but with a modified potential 

V {v) {x;a) = V(x;a) + 2w' v {x;a). (5) 

This can be schematically resumed as 



W\ >-> w\' 

A(w v ) 

V >-> v( v h 



(6) 



The corresponding eigenstate of H^ v \a) = —d 2 /dx 2 + V^ v '{x; a) can be written 

. cxp (- /^r,„,) ~ ' A ( ^ M , P) 

V J J V Ex (a) - E v [a) 

where A (a) is a first order operator given by 

A (w v ) = d/dx + tii„(i; a). (8) 

Eq([7]) can still be written as 
where W (y\, j/2 | x) is the wronskian of the functions y%, 1/2 ■ 
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From V, the DBT generates a new potential V^ u ' (quasi) isospectral to the original one and its eigenfunctions are 
directly obtained from those of V via Eq©. Nevertheless, in general, w u (x;a) and then the transformed potential 
are singular at the nodes of ip u (x;a). For instance, if tjjn(x;a) [y = n) is a bound state of H(a), is 
regular only when n = 0, that is when ip n —o is the ground state of H, and we recover the usual SUSY partnership in 
quantum mechanics [22l . [23| . 

We can however envisage to use any other regular solution of Eq© as long as it has no zero on the considered real 
interval /, even if it does not correspond to a physical state. For some systems, it is possible to obtain such solutions 
by using specific discrete symmetries Ti which are covariance transformations for the considered family of potentials 

{ r a ^ ai (10) 

[ V(x; a) 4 V(x; a t ) = V(x; a) + U (a) . 

Ti acts on the parameters of the potential and transforms the RS function of a physical excited eigenstate w n 
into a unphysical RS function v ni i(x;a) = Ti (w n (x;a)) = w n (x;ai) associated to the negative eigenvalue £ n j(a) = 
Ti(E n (af) = U (a)-E n (ai) < 0. ' 

- v ' n ,i( x ; a) + v l,i( x > a ) = v ( x 't a ) - £n,i( a )- ( U ) 
To v n> i corresponds an unphysical eigenfunction of H(a) 



<p n ,i(x; a) = exp I - J dxv n ,i(x; a) J (12) 

associated to the eigenvalue E n> i{a). 

If the transformed RS function v ni (x;a) of Eq (|lip is regular on /, it can be used to build a regular extended 
potential (see Eq© and Eq©) 

V^{x;a) = V(x;a)+2v' nii (x;a) (13) 
(quasi)isospectral to V(x; a). The eigenstates of V^ 1 ' 1 ' are given by (see EqQ) 

w { k n ^(x;a) = -v n ,i(x;a) + v ^\i%^Xi) (14) 
^■ l \x;a) = exp f dxw ( ^' l) (x; a)) ~ -j=J= = A(v n ,i) ip k {x; a), 

for the respective energies Ek(a). 

The nature of the isospectrality depends if 1 / <p n ^(x; a) satisfies or not the appropriate boundary conditions. If it 
is the case, then 1 / (f> nji (x; a) is a physical eigenstate of H^ n ' l ^{a) = —d 2 /dx 2 + V^ n,l '(x;a) for the eigenvalue £ n ,i{ a ) 
and we only have quasi-isospectrality between V(x; a) and P^"' l )(x; a). If it is not the case, the isospectrality between 
V (n ^(x;a) and V(x;a) is strict. 

The above construction can be summarized by the following diagram 

Ek > : Wk >— ^ w i n '^ }Physical RS functions 

£k,j < : Vk,j ^ ' v k!j^ } Regularized (unphysical) RS functions 



with 



V >-> V (n '^ Potentials 



= A(v n ,i) ip k 

&k\j = A(v ntl )(f> kd . 



(16) 



This procedure can be viewed as a "generalized SUSY QM partnership" where the DBT can be based on excited 
states RS functions properly regularized by the symmetry Tj. 
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III. DISCONJUGACY AND REGULAR EXTENSIONS 

To control the regularity of v n ,i we can make use of the disconjugacy properties of the Schrodinger equation for 
negative eigenvalues. 

A second order differential equation like EqflTJ is said to be discord ugated on 7" C R (V(x; a) is supposed to be 
continuous on I) if every solution of this equation has at most one zero on I [U [28| . As it is well known, this zero 
is necessarily simple and at this value the considered solution changes its sign. For a closed or open interval I, the 
disconjugacy of EqJT]) is equivalent to the existence of solutions of this equation which are everywhere non zero on / 
[I?], HH. In the following we will consider I = ]0, +00 [. 

We have also the following result: 

Theorem [U HH If there exists a continuously differentiable solution on I of the Riccati inequation 

- w'(x) + w 2 {x) +G(x) < 0, (17) 

then the equation 

^"(x) + G(x)i>(x) =0 (18) 

is disconjugated on /. 

In our case, since £ n ,i(a) < 0, we have 

- w' (x; a) + w 2 (x; a) = V(x; a) < V(x; a) - £ n ,i{a), (19) 

Wq(x; a) being continuously differentiable on I. The above theorem ensure the existence of nodeless solutions <p(x; a) of 
EqlJU) with E\(a) = £ n ^{a), that is, of regular RS functions v(x; a) solutions of Eq([TTj). To prove that a given solution 
4>(x;a) belongs to this category, it is sufficient to determine the signs of the limit values 4>(0 + ;a) and 0(+oo;a). If 
they are identical then <f> is nodeless and if they are opposite, then <f> presents a unique zero on I. 

In the first case V (x; a) + 2v'(x; a) constitutes a regular (quasi)isospectral extension of V(x; a). 

In fact, Bocher [23] has even established a more precise theorem: 

Bocher 's disconjugacy theorem 

If we can find u continuously differentiable on [a, b] such that 

u'(x) +u 2 (x) < G{x) (20) 

on all [a, b] , then nor ip(x) — u(x)ip(x) — ip'{x) nor ip(x) can have more than one zero in [a, b]. Moreover, if ip(x) 
vanish in one point of [a, b] then tJj(x) is everywhere nonzero in this interval (and reciprocally). 

If this disconjugacy property is satisfied for every [a, b] c ]c, d[, then it is in fact satisfied in all the open interval 
}c,d{. 

Consider first the case where we have G(x) > 0, Va; € ]c, d[. We can then take u(x) = which implies that any 
solution ip{x) of Eq (|18[) as well as its derivative ip'{x) have at most one zero in ]c,d[. If £0 is the corresponding zero 
of ip'(x) then ip(x) is strictly monotonous on ]c, and ]^o, d[. 

If we consider now the case of Eq([T]) on I = ]0, +oo[ for a strictly negative eigenvalue £ n ^{a) < 0. From Eq ([T^|) we 
can take u{x) = —wo{x;a) to apply Bocher's theorem from which we deduce that not only Eq([T]) is disconjugate on 
]0, +oo[ but also the SUSY partner V^ ' (x; a) of any solution of Eq([T]) can have at most one zero on ]0, +00 [. Indeed 
(cfEq©) 



(21) 
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IV. ISOTONIC OSCILLATOR 

Consider the isotonic oscillator potential (ie the radial effective potential for a three dimensional isotropic harmonic 
oscillator with zero ground-state energy) 

V (x- u, a) = ^-x 2 + a(a ~ 1} + V (u, a), x > 0, (22) 
4 x z 

with a = I + 1 > 1 and Vq(ui, a) = —u> (a + ¥)■ It is the unique exceptional primary translationally shape invariant 
potential of the second category [2l[ . The corresponding Schrodinger equation is the Liouville form of the Laguerre 
equation on the positive half-line and its physical spectrum, associated to the asymptotic Dirichlet boundary conditions 

ip (0 + ;u,a) = = ip(+oo;oj,a) (23) 

is given by (z = uix 2 /2, a — a — 1/2) 

E n (a) = 2nu, V„ [x; to, a) = x a e- z ' 2 L« (z) . (24) 
To ip n corresponds the RS function [2l[ 

w n {x; u>, a) = Wo{x; u>, a) + R n (x; u>, a), (25) 

with 



and 



I \ U a ,^„s 

w {x;uj,a) — — x (26) 

2 x 



, \ -2nu , , 2 On — j + 1) u , . 2uj 



ujx — (2a + 1) /x— <jJX — (2 (a + j) — 1) fx— lux — (2 (a + n) — 1) jx 

uxL a n +_\{z)/L a n {z). (27) 



The shape invariance property satisfied by V is [2lM24| 

V(x; w, a) = V(x; u, a) + 2w' (x; u, a) = V(x; w, a + 1) + 2oj (28) 

We have three possible discrete symmetries for V (see Eq (fT0|) ) which are given by 
1) 

u r 4 (- W ) , { V(X] U l °) ^ V{X; W ' a) + U{2a + 1} (29) 



2) 



3) 



w n (x;uj,a) -4 v ni i(a;;a;,o) = -w, a), 



a £_ 1 _ a f y (a! ; W| «) ^ l/(x; W| o) + a,(2o - 1) (3Q) 
{ w n (x;oj,a) -> u„ )2 (»;w,a) = w n (x;u), 1 - a), 



(a,, a) r3=I 4° r - (-a,, 1 - a) { V ^ W ' a > ^ V ^ W ' a ) + ^ (31) 

( w„(x;w,a) -4 v ni3 (x- 7 uj,a) = w n (x; —u, 1 - a). 
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In the {oj, a) parameters plane, r+ and T_ correspond respectively to the reflections with respect to the axes oj = 
and a — 0. The RS functions v„.i, i — +, —, 3, satisfy the respective RS equations 

- v 'n,i( X > w > a ) + V l,i( X > w > a ) = V ( X ' w > a ) ~ a ) ) ( 32 ) 

with 

£ n>+ {oj,a) = £J_( n+a+1 /2)(w) <0 

£ n ,-{w, a) = £ n+1/2 _ (w) (33) 
£ n , 3 (a;,a) = E_(n +1 -)(w) < 0. 

These eigenvalues are always negative in the i = + and i = 3 cases and the inequality £„ r —{oj, a) < necessitates 
to satisfy the constraint a = a — 1/2 > n. When £ n ^{oj, a) < the Schrodinger equation 

<fi"(x; oj, a) + (£ ni (oj, a) — V{x; oj, a)) 4>{x; oj, a) = (34) 

is disconjugated on ]0, +co[. We are then sure that everywhere non zero solutions of Eq (l34p exist, the question being 
to if <fi n .i(x; oj, a) — exp (— J dxv n .i(x; oj, a)) is such a solution. 

V. LI AND L2 SERIES OF EXTENSIONS AND CORRESPONDING EXCEPTIONAL LAGUERRE 

POLYNOMIALS 

The LI and L2 series of extensions is obtained using the T + and T_ symmetries respectively in which case we have 
synthetically 

£ n .i {oj,a) = -2oj (a + i (n + 1/2)) s = ± (35) 
tf>n,i(x;u),a) = ^o,-i(aj;w,a)i/^ a {-iz) ' 

with 

4>o,i{x; oj, a) ~ x la+1/2 exp (iz/2) , (36) 
where a is supposed to statisfy the constraint a > —in. The corresponding RS function is given by 

v n j(x; oj, a) = v ,i(x; oj, a) + Q n j{x; oj, a), (37) 



with 



and 



1 \ - u ta + 1/2 1 ... 

v .i(x; oj, a) = -i—x = (1 (z + a) + 1/2) (38) 

2 x x 



Q n ,i(x; oj, a) = -iuxl££ {-iz) /L™ {-iz) . (39) 



Since 30, 31 



[ft /'™\ (~1) ~ra 

a:— J-+oo 



(40) 



where (X) = (X) ... {X + n — 1) is the usual Pochhammcr symbol [3l|, we have for a > 1/2 



(+00; w, a) = +00 
t> n ,+ (0+;u;,a) = 0+, 



(41) 
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and for a > n 

4>, h+ (+00; to, a) = ± 



+ (0+;w,a) = ±oc, ( 42 ) 
with ± = (-1)™. 

Because of the constraint a > —in, Eg ([34]) is disconjugated both for i = + and i = — . Combined to Eq flTj) 
and Eq (|42[) . this ensures that (f>n.i keeps a constant strictly positive sign on ]0, +oo[ which means that v ni (x;uj,a) is 
regular on this interval. Note that this result implies in particular that is L" (x) nodeless on the negative half line and 
that L~ a (x) has no zero on the positive half line when a > n, which is in agreement with the Kienast-Lawton-Hahn 
theorem (30L l3lj|. 

V( n '*)(x; u, a) — V(x; oj, a) + 2v' n ^x; u), a) is then a regular extension of V(x; w, a). Since l/cf> n ,i (x; u), a) diverges 
at the origin, it cannot be a physical eigenstate of H^- n ^ . Consequently, V^ n ^ and V are strictly isospectral. 
The (unnormalized) physical eigenstates of V^- n ' l \ ^™'^(o;;u;,a) — exp f— J dxw k n ''^ (x; U), a)^j , satisfy 



(w, a)^ l) (x; uj,a) = E k («) ^ A (x;oj,a) 
u { k n ' l) '(x;uj,a) + (w { k n ^(x;u;,a)y = V^\x;uj,a) - E k («) 



(43) 



where, using the shape invariance property of V (cf Eq (j2"8)l ). we can write 

F (n ' j) (a;;w,a) = V(x; u, o_i) + 2Q' ni {x; u, a). (44) 

From Eq (rj"4")) we have also 

4> k n,l \x;u),a) = A(v n;i ) ip k (x;u,a), k>0. (45) 
More precisely, using the identities [3(J HU 



LI (x) + I%±\ (x) = L- +1 (x) 
xL^+l (x) - aL«_ x (x) - nL^ 1 (x) = 0, 



(46) 



we obtain 



-0^™' l) (a;;w,a) = (v nii (x; ui, a) - Wk[x; w, a)) ipk(x; u, a) (47) 
i x a +* exp(-z/2) 

where 

Kk,a (z) = (k + n + a) LI (z) L~ a (z) - {-n + a) L% (z) L~^ (z) - (k + a) L^_ x (z) L~ a (z) , 



(48) 



are polynomials of degree n + fc, namely exceptional Laguerre polynomials (ELP) of the LI and L2 series respectively. 
From the orthogonality conditions on the eigenstates for fixed values of n and a, we retrieve the fact that the ELP 
L l nka (z) constitute orthogonal families with the corresponding weights 

"«(«>- ,49) 
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VI. FORMAL SCHEME FOR MULTI-STEP DBT 



A. General elements 



The question of successive iterations of DBT is very natural and is at the center of the construction of the hierarchy 
of hamiltonians in the usual SUSY QM scheme [4l[ . Staying at the formal level, it can be simply described by the 
following straightforward generalization of EqQ 



(50) 



where yy£ Ut "* Urn ' is a RS function associated to the eigenvalue E\ of the potential 



m — 1 



I/l„...,l^_l) 



(x;a) 



3'=1 



The corresponding eigenfunction is given by (cf Eq(0 and Eq®) 



(x; a) = 1 ^ 1 --^- 1 \x; a) = A (wfc-">^) ..A^fafaa), 



that is, 



w 
W 



(4. 



where TU (yi, y m | x) is the wronskian of the functions y±, y„ 



W(yi,-,2/m I a;) 



?/i (x) ... y m (x) 

(m— 1) / \ (to— 1) / \ 

Vl (X) ... ym (x) 



(51) 



(52) 



(53) 



(54) 



Other useful representations for the extended potentials and of their cigenfunctions are possible. From Sylvester's 
theorem we can deduce [361 ] the following property for the Wronskians 



-,y m ,y | x) 



W (W (yi, ...,y m ) ,W (y!, ...,y m -i,y) I x) 
W (yi,...,y m -i | x) 



(55) 



Using the identity [3C 



W(uyi, ...,uy m | x) = u m W (yi, y m | x) , 



(56) 



this gives 



W (?/!,..., y m , J/ | x) _ ^W(vi,...,y ro _i)' 



,y m -i,y) 



,Vm-l) 



W(yi,...,y m | x) 



VV r (2/i,...,!/mb) 
W(yi,...,y m -i\x) 



(57) 



and comparing to Eq (|53p . we obtain 
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I (v 1 ,...,v m ) f \ 

[X ' a) - W^ Vl ,...,^ Vm \x) ' 

which is the well known Crum formula (37l . |40| for the eigenfunctions. 

Inserting this result in Eq fSTj) . we then deduce the Crum formula for the potential [37l lie 



(58) 



V ^-^(x;a) = V(x;a)+2(logW(iP„ 1 ,...,ip Um \ x))" . (59) 

Since all the functions ij) v . implied in the wronskians in Eq (l58l) and Eq (|59[) are eigenfunctions of the same hamil- 
tonian H(a), the properties of the determinants allow us to replace in these wronskians the even derivatives ipffl by 
(—E Vj ) ip Vj and the odd derivatives i/>£f +1) by (-E Vj ) 



For instance 



where 



and 



where 



W (i/j Vl ,...,ip U2l | x) 



[x; a) 
(x; a) 



(-E Vl ) 1 ip„ 1 (x;a) 



ip„ 2l (x; a) 
i>l 2l 0; a) 



21 



(x; a) 



\Z-l 



(E Vl ) w Ul (x]a) 



tj) vx (x; a) 
ip' Ul (x; a) 



(-E Ul ) 1 ij)' Vi (x;a) 



w U2l (x; a) 

(E.J- 1 
(E V2l ) 1 1 w U2l (x;a) 



ij) V2l (x] a) 
f V2l (x] a) 



(-Eu 2 i +1 ) ip' U2l+1 (x;a) 



{-E Vl ) ip Vl (x;a) . 

21+1 

= A^>-> v *>+^(x;a) l[ i> Vj (x;a) 



(60) 



(61) 



(62) 



1 

w Vl (x; a) 



w V2l+1 {x; a) 



(E V1 ) w Ul (x;a) ... (E Un ) w U2l+1 (x; a) 

(E Vl ) ... (Ey 2l+1 ) 



(63) 



Then we obtain the representation of the ^" 1 ''"' l ' m ' ' m terms the so-called Crum-Krein determinants ^i v u—i^m) 

M M M 



<P X (x;a)=i> x (x;a) — yj—^ . 



(64) 
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Until now, all these results are purely formal. The central problem is now to choose in an appropriate way the 
family of eigenfunctions (VVu •••) VVm) from which the successive DBT are built (see Eq ([55|) and Eq (|5Tj) ) in order to 
ensure the regularity of the successive extended potentials yC" 1 '— > I/m ). 

A first answer has been given more than half a century ago by Crum [37| . The Crum proposal [37J is to take the 
successive eigenstates of the discrete bound states spectrum starting from the ground state 

(n, ...,u m ) = (0, ...,m-l) (65) 

Krein [H, HU, EU extended this result showing it is possible to choose sets of unnecessarily successive eigenstates if 
they satisfy a given condition, namely 

(v u ...,v m ) = (m,...,n m ) e N"\ (66) 

with 

(n — ni){n — 712). ..(n — n. m ) > 0, Vn G N. (67) 

This last is verified if the set (n 1; n m ) is constituted by "aggregates" of an even number of eigenstates. The 
possibility to employ sets of two juxtaposed eigenstates has been rediscovered later by Adler |39| . These results have 
been used extensively in the context of higher order SUSY (see for instance [25L l4l| - l52| ). In (48[, Samsonov has 
still extended the Krein- Adler result. Note finally that Fernandez et al (4|| employed the specific "Wick rotated" 
eigenfunctions of negative energies introduced by Shnol' [IH 13 E3| to build successive extensions of the harmonic 
potential. 

B. Application to the isotonic oscillator 

Our aim is to consider the possibility to generalize our construction in a multi-step version on the basis of the 
regularized RS functions (or unphysical eigenfunctions) obtained above. Staying at a formal level, it would correspond 
to the following generalization of the diagram EqlfTS"]) (see Eq(|50|)) 

A(«„, u) , ■ , M^%°) , ■ ■ » A( „(»l.<l" ".»m-l.<m-l) ) 

E k >0: w k ~ w k nun) ~ w (niM,n 2 ^) " ^ |u (»i,ti,,...,W B 

ru 

A(v ni n ) . . , A(„££>) A(v (ni,^.,...,n m _ lttm ^) ( 68 ) 

c n v (ni.ti) 2, 2 (ni,«i,n 2 ,t2) ™m..«m (n 1 ,% 1 ,,...,n m ,i m ) 



A(t) («l.n..-.»m-l.<m-l) ;) 



y(ni,ii„...,n m ,i m ) 



with ?4 = Wfe. Under this way, we generate chains of isospectral extensions V^ ni ' il "'"' nm '' l " m ^ of the potential V with 
(see EqdSU and Eq([59]l) 

F (nx,u,...,n TO ,i m ) (;z , ;a) = y ^. a) + 2 ^ 1) n,..,n m _ 1)im _ 1 ) (a;;a) y (69) 

m 

= V(x; a) +2 (log W ..,<A nm ,i m | x))" 

w (ji 1 ,i 1 ,,...,n m ,i m ) ^ e fni,n„...,n m ,i m ) j^s functions of the extended potential y( n i' 11 '— '"""'•») associated 
respectively to the eigenvalues Sfc and £ k ,i- 

The corresponding eigenfunctions are given by (see Eq([53| and Eq(f58jl) 
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} (x;a) 



„,/ "m-li^m-l) , (»li<l "m-lilm-i; i \ 



W^K.n.-^m-l.im-l) ,(»l.<l.....«m-l.im-0, ^ [ ' ^ > 

(ni,ii,...,n m> i ra )/ . x ^ |x J _ ff^.n ^, im ,» t|i |») 



Suppose that at the step to — 1 all the extended potentials in the chain are regular and strictly isospectral. The prob- 
lem is now to control the regularity of the potential V^ ni ' n '"' ,nm ' lm ^ and of its physical eigenstates ^ ni > ll "---> nm - 1 '* m - 1 ) 
as well as its strict isospectrality with y(™i> l i,---,nm.-i,4 m -i)_ This is achieved if the unphysical eigenfunction 
& { Ztj''' nm ~ 1 ' im ~ 1 \ associated to the DBT A(«^^"'"' nm_1 ' im "^), is nodeless and if i/^.£»-". n »-i.*m-0 can " 
not be into the set of physical eigenstates of V^ n ' 1 ' %1 '"'' nm ~ 1 ' tm ~ 1 ' . 

The disconjugacy of the Schrodinger equation 

f{x;u,a) + (£„ m ,i m (u,a) - V ( n ^-^-^-^(x;oj,afj cj>{x;u,a) = 0, (71) 

for a negative eigenvalue £n m ,i m guarantees that w^'^''"'"" 1-1 '*"* -1 ' is at most singular in one point. The regularity 
of the extended potential ^("i.'i^'^.'J and its strict isospectrality with the preceding one in the chain are then 
satisfied as soon as the unphysical eigenfunction (fr^'i 1 " "' , ™ m ™ 1,lm ~ 1 ^ satisfies appropriate boundary conditions, which 
was already the argument used above in the one step case. Namely, the strict isospectrality and the regularity are 
satisfied when 0(" 1 ^ 1 '--- ,nm - 1,lm - 1 ) tends to oo at one extremity of ]0, +oo[ and to at the other extremity, with the 
same sign for both limits. 

In the case of the isotonic oscillator, there are in fact three type of regularization transformations £ = +,—, 3 
associated respectively to the three series of EOP LI, L2 and L3 [li|. At each step we can a priori use any of the 
three type if regularized RS function. However, since the use of DBT associated to RS functions of the L3 series leads 
only to quasi isospectrality, the argument above is no more valid. Consequently, in the following we exclude this case 
and consider only DBT built from the RS functions of the LI or L2 series. 

From Eq ([64| . we see that in the present case 



0£±"" , ' ,m " l,<m ~ l) tew.a) = <Po,i m (i;^)^""' 1 "™'' J fow.a), (72) 

where 

A(m,ii„...,n m ,i m ) / \ 

R {m M „.. Mm ) ( a) = u a ( _ • } . 7'"'^ - (73) 

v ; n m V m I ^(ni ) ti„... ) n m _i ) t m _i) ( X ;uJ,a) 

is a rational function. Due to the presence of the exponential in 0o,i m , the behaviour at infinity of <^™ 1 '* 1 '---' rim - 1 < lm - 1 ) 
is the same as for 4>o,i m , namely 



^(niA,...,n m _ ui m - 0^.^^ „ exp(z/2)x a+1 / 2 R( n i^»-> n ™-+) (+oo;uj,a) -> ±oo 

m ' x— J-+00 X— >+O0 

^ ni ^'-' nm - l ' <m - l) (a:;a;,a) ~ ex P (-«/2)x- 0,+1 /aij(«i,ii,,-, "-.-) (+oo;w,a) -> 0±, 

m 1 2;— »+og >-+oo 



(74) 



(the sign of the limit being the one of RS ni * 11 »••■>"<».+) (+oo;w,a) or .R(™i> 4 i"---> n ™>-) (+oo;cj,a) respectively). 

The absence of node for 'I 1 '"' ,nm - 1,lm_1 ^ and the strict isospectrality between yffi^.-^m-i.im-i) a nd 
^(ni,ii,...,n m ,i m ) ) are simultaneously ensured if and only if we have 

</»t!;+ 1 *-* nm - 1 * im - l) (x;a;,a) ~ a .a+ 1 /2 fl (n 1 ,i 1 „... J n m ,+) (0+; W ,o) ^ + 0± 
^(m.n.....n™- l ,*m-i) ( j 1"* - B -«+l/2jj(ni.ii„...,n m ,-)(0+; W ,o) ^ ±00, ' ' " ' 

(the sign of the limit being the one of i?( n i>*i»— >""»,+) (o+;aj,a) or iZ^ 1 *' 1 """ >™ m > _ ) (0 + ;w,a) respectively), with 



si<?n ^(m.*i..-.n»,im) (0+;w,o)) = «ffti f Jj(»i.*i»-.«»A») ( +0 o;u,,a)) 
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VII. TWO-STEP CASE 

Before to consider the general case, we first consider the two-step case. We then have to check that the unphysical 
eigenstates (f>^^ satisfy the appropriate boundary conditions. Using Eg ([55)) in Eg ([7D|) . we have for m = 2 



An u ii)( \ _ W ( X ' W ' a ) L ni ; (gi ^> «)^n 2 " (-^z) \ x) 

< Pn 2 ,i 2 ( X ' UJ ' a ) — , / \Tiiat ■ ^ 



where a > —ijUj, j = 1, 2, or with Eq (l56p 

(z; w, a) = o ,u (a;; w, a) 

If ii = %2 = i, then 



W {-hz) , (x a exp (z/2)) 12 41 (-i a *) | x) 

Lli" [-hz) 



(76) 



Using [36 



and 130, 31 



this gives 



W(LZ(-iz),L™(-iz)\x) 



( x > w > a ) = ^o,i (»; w, a) " V ^ ni - YT^riV\ 



TU(yi,...,y m | .t) = ( — J W(y x ,...,y m \ z) 



(LZ (x))' = -L%±l (x), 



( i} iwx* a+3 / 2 exp(iz/2) 

w to u > a ) = 



LZ i-iz) 



L% (-iz) L% (-iz) 



Lzt\ Hz) v::a {-iz) 



n 2 
za+1 



We then have (cf EqflDJ)) 



4& 



<^Si l) (^w»a) *" 2 («2 - ni) 



/n 1 +ia\ /n 2 +ia\ 
\ rti — 1 / V n2 — 1 / 

_^ x 2n 2 +fa-l/2 exp( -^ 2 ), 



V n 2 — 1 / n 2 



that is, 



(77) 



(78) 



(79) 



(80) 



(81) 



V ™2,+ V , i ) x ^ 0+ \ n2 -l) n 2 x ^ Q+ 



x—>-\-oo 



where ± = sign(ni — m), and 



±oo, 



if a > 1, 



x— f+oo 



(82) 



/(ni, 



} (a;:w,o) - (-if 2 W x-(«- 2 )-V2(«-" f ) 2a=2L 



n 2 , 



£_ _) (a;;w,a) (-1)™ 2 f^-^ll ^-a-1/2 exp ( _ z/2) 

where ± = ( — l)™ 2 sign(ri2 — n±). 



-> ±oo 

' "" , if a > sup (ri2, 7ii) 



(83) 
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If i\ = —12 = —i, then using the last of the following identities [5J 



(z a L« (z))' = (n + a)z<*- 1 L%- 1 (z) 

{e-*I% (*))' = -e- z ^ +1 (z) (84) 
(z a e-*L% (z))' = (n + 1) s^e-I^ 1 (z) , 



we obtain 



'„ 2> i (^;w,a) = m +/ exp (-«/2) I - J ■ 

L-™ (iz) zL% (-iz) 



(85) 



2x ia-l/2 exp (j z / 2 ) 



We then have (cf Eqffl) 



^^fow.a) ~ i "^(|)" 2+1 ^+»+ 3 / 2 exp(zz/2). 



(86) 



Consequently (a > 1) 



(87) 



and when a > 



with ± = (-1)™ 2+1 . 

When a > sup(ri2, ni), we have the correct asymptotic behaviour which ensures the regularity of V {ni ' il ' n2 ^(x; cj, a) 
and V^ ni - ll \x;uj,a) as their strict isospectrality with V(x;ui,a) for every choice of (ni,ti, 712^2) with ni 7^ ri2- 
We then obtain a chain with two strictly isospectral regular successive extensions of V 

A(v„ i ) ^("i" 1 ^ 1 ') 
y 1 y(m,n) y(rn,i I ,n 2 ,i2)_ fgg) 

VIII. GENERAL M-STEP CASE 



Consider now the general case of a m-step DBT. We proceed by induction, supposing that we have built a chain of 
to — 1 regular strictly isospectral potentials 

A( Vn , ) ^i" i r i> ) A(4" i : < 1 i T , : n r" 2 ' im_2) ) 

y >_V 1 y(ni,n) ^ 2 ^ m ' m M y(ni,«i,...,n ro _i,j m _i)_ j-qq-j 

and that we have a > —ijiij, Vj < to. This implies in particular that the cj^™ 1 '* 1 '''''™* -1 ' 1 * [x\uj,a), J < to — 1, 

which are eigenfunctions of if ( n i>*i>— > n, »»-ii*m- 1) associated to the respective negative eigenvalues £ rlj _ 1 ,i i _ 1 {uj, a) < 0, 
are all nodeless on ]0, +oo[ and that they keep a constant sign on this interval that we can always take as positive. 
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We have to determine what is the constraint to implement for the unphysical eigenfunction ^C™'* 1 '•••'"""> 1 ' lm x ) m 

order that it satisfies the good boundary conditions necessary to build a regular and strictly isospectral extension 
y(n u i u ...,n m ,i m ) at the next step _ 

From Eq ([74|) we know that 



, (ni,ii,...,n m _2,i m -2) / \ 

b n m U+ J {x;w,a) -> +00 

x— >-\-oo 



(91) 



and consequently we have also (cf Eq (17510 



% m li+ '(aw, a) 

, (ni,ii ,...,n m -j,t m _2) / \ 

r n ' _ '{x;u,a) -> +00. 



(92) 



For x sufficiently large or in the vinicity of 0, we have £ nrn _ lt ±(u),a) < V^ ni, ' ll '"'' nm ~ 2 ' lrn ~ 2 \x; oj, a). The Bocher 
disconjugacy theorem then implies that in the neighbourhoods of and of +00, the functions < ^ 1 ^<^ , <™™-*< im - 3 > are 
strictly monotonous, that is strictly decreasing if i m -i — ~ and strictly increasing if i m -\ = +• 

We already know (cf Eg ([73])) that 

, (ni,ii,...,n m _i,i m _i) / \ J ioo, if i m = + /r\o\ 

%i m ,i m '{+oo;w,a) = i o ± , ifi m = - ^ ' 

We prove now by induction the following result for the unphysical eigenfunctions of the potentials obtained at the 
(m — l) th step (z = ujx 2 /2) 



^— nm - ulm - l] {x-uj,a) ~ ^^'•••^-i' i --0( y m (a+ g (ii,..,i m ))-l/2 (1 + o ^ ))) (94) 

m x _ TTl, TT1 

is the same sign as <p£ 1 '] 1 ''"' nm_:L ' lm ~ I '(+oo;a;, a) and where q(ii, ...,i m ) = 
q + (ii, i m ) — q_ (ii, i m ), q± (ij, i m ) is the number of ij equal to ± in the set i m ). g± i m ) can be 

viewed as the number of state of "± charge" used in the chain of DBT at the m th step and the difference q (ii, i m ) 
as the total charge associated to the chain. 
We then have 

z (ni,2i ,.. . ,n m _ 1 ,i m — 1 ) / \ r\4- -r r • ■ \ , -i / ■ \ 

^ , '{x;w,a) -> 0^, if a > -g («!,..., + 1 = -g z m _i) 

, (ni,ii,...,n m _i,i m _i) / n . • r / . ■ \ , / ■ s \ aKj J 

r nm - '{x;u,a) -> ±00, if a > -g (ii, z m _i, -) = 1 - q (ri, i m -i) , 

which are necessary conditions in order that the potentials y( n i,Ji,---,« m ,+) an d y{n 1 ,i 1 ,...,n m ,-) reS p ec ti V ely are regular 
and strictly isospectral to V^ ni,%1 "'" , " m " I,,m " 1 ', Note that these constraints are independent of the rij but depend 
only on the type of the states used to build the chain. 

Refering to the results obtained in the 2-step case (see Eq ([82|) . Eq ([83| and Eq(|88j)). we can verify that this property 
is valid for m = 2 since 

9 (+,+) = -?(--) = 2 (m) 
Suppose that Eq (IM|) is verified till the (m — l) th step, that is, 



$± il '-"' nm - 2 ' im - 2) {x;u,a) ~ a;±(a+ 9 (ii,..., *™-2»+i/2 (1 + 0(2)), (97) 

with a > -g(ii, ...,i m _2, +) = -1 - <? («i, ...,i m _ 2 ), if i m _i = +, and a > 1 — q (i l5 ...,i m _2, -) = 2 - ...,i m _ 2 ), 
if 4-i = -• 
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Since 2 ' %m ^ anc ^ 0^" 1 '* 1 ''"' n * B 2 ' 1 '™ 2 ' are both eigenfunctions of H^ ni ' n, --- ,nm - 2 ^ m - 2 \ using the wron- 

skian theorem [55l | in Eq (|70[) we can write the unphysical eigenfunctions at the (m — l) th step as 



i;w,fl = 7 — : : ; (-98) 

TJ7 /" i(ni,il,».,n»-!,ii»-!) j(Tll,il,...,7lm-2|tm-2) I \ 

W \fn m U,i m ^ >0fe,i 1^0 J 



,(m,ii,...,n m _i,i m _i) 
n m -l>*m-l 



(x;u},a) 



We have to consider several cases corresponding to the four possible couples (i m _i,i). 

1) Consider first the case i = i m -x = + in which we take xq = 0. Using Eq ([92|) . Eq ([98]) gives 



r x Ani,ii,...,n m -3,i m -a) ft \ Ani,ii,...,n m -2,i m -2) ft \ 

.(n 1 ,i l ,...,n m _ a ,i m _„n m _ 1) +), * Jo ^n m -i,+ (&^a)<ftfc,+ g^Mp 

fflL, (xtu;, a)~ -, — : : > . (99) 

^n m _i,+ (a; a;, a) 

Since and ^H,..,nm- 2 ,i m -2) are both gtrictly positiye Qn ]o 7+00 [, ^A..."....^^.^ 

keeps a constant sign on this interval and we can always take it as positive. Using Eq (l97l) . Eq ([99|) gives 



™'"- 2, '"- 2 '"'"- 1 ' +) (x;^a) ~ x (a+ q (i u ...,i m -2,+,+))-i/2 (1 + ? (100) 

Consequently, in this case, Eq dMl) is also verified at the m th step. 

2) In the case i = i OT -i = — , if we take xq = +oo, we have similarly with Eq fliTTI) and Eq (j9"5)) 

p+oo , t Am, ii,...,n m - 2, irn -2) ft v i(ni ) «l,... l n„_j,» m _ s )/ t x 

<pl '(x;u,a)~ -. : : r . (101) 



As in the preceding case, since 4>n^l\'.-' nm 2,lm 2 and ^j^'* 1 '— >" m 2 ' 2m ^ are both strictly positive on ]0,+oo[, 



<,,... , _ ^,iu-,n m -2,i m -2) anA Jnuiu-,n m -2,i m -2) 

^ni,ii,...,n m _2,i m -2,n m _i, ) k ee p S a cons tant sign on this interval and we can always take it as positive. From Eq(|74|) 
we deduce that for any e > 0, dfaj>!^ 1 '*l''l'' nm ~ 2 ' lm (fa oj, a )0(™^ , * 1, "-> n " l - 2,l ' n - 2 ) (fa W; a ) j s finite and positive. 
Since (a > 1 - g (ii, i m _ 2 )) 



j (ni,ii,...,n m _2,i m -a)/„. , , ,iii- 1 »m-2 I im-j) 
^n m — 1,— 



(x; w, a^^-^-2^-2) ( x . Wj a) _ a .-2( Q+9 (i 1 ,...,i m _2))+i (1 + 



Eq fTUTI) gives 

0(ni.ii.....n m - a ,< m - a ,n m _ 1 ,-)/ j _ a .-(a+ 9 (u ) ... ) i m - a ,- ) -))-l/2 - 

Eq (|94l) is again verified at the m th step. 

3) It remains to consider the cases i = + = —i m —\ and i = — = — i m -i- Since 0^_' ll ''"' nm - 2,lm - 2 J an( j 
^(ni^n,...,n m -2,im-2) nave n0 ^ ^he same behavior nor at nor at +00, Eq dMl) is not adapted and we proceed to a 
direct evaluation of the wronskian W ($^'2l' l -' nm ~ 3 ' tm ~*\ <^™^ , * 1, '" ,nm_2, * m ~ 2 ) | xj. 

Since ^(^> l i>--->" m - 2 >' i "'- 2 ) anc i ^(ni.»i.'-.nm-2,« m -2) arg s ^ r j c ^iy positive and respectively decreasing and increasing 
both at infinity and in the vinicity of the origin, W ^^^ L 'l'' n "*~ 2 ' lm_2 \ 0^™|' n '---' nm - 2 ' i ™- 2 ) | x ^ the same positive 
sign at these two limits. Due to the disconjugacy of 

<p"(x;u,a) + (e ki+ (u,a) - V*" 1 '* 1 --""*- 1 '*"-^!; u, a)) (f>(x;u,a) = 0, (102) 
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we then deduce that ( p(^ ll >—> nm 2 > lm 2 >™ m 1 > }(x;u),a) keeps a constant sign (taken as positive) on all the interval 
]0,+oo[. 

Using (see Eq ([5?]t ) 



/ (ni ,ii,...,ri m — 2,i, 



(a:;a;,a) ~ a; («+9(n,..^-2))+i/2 n — 2 ^- 2) (0) + O(z)) 

(s;w,a) - i-(a+9(<iv..,U- a ))+i/2^i.n. r .«m-».*™-a)(o) + 0(« 



(103) 



we obtain 



j^r ^(m.U,— ,n m _2,i m _2) ^(wi,*i,-..,n m -2,»m-a) | 

~ 2(a + q (t x , z m ^ 2 )) fc«-»,»»-.,<»-,) (0 )giA "»-»,<»-.)( ) + ( 2 ). (104) 

a:— >0+ 



Consequently, we deduce 



and 



6 (ni,ii,...,n m _ 2 ,i m _ a , nm _ 1 ,-) (a . x _ a .(<H-9(<i,...,<m- a ,-,+))-l/2 (105) 



n/fltnii'i.-inm-Jiim-j) ,(rn,ii r ..,ti m -2,i»,-2) i \ 
A- \x ] u: 1 a) = ,(„ 1 , n ,...,» m _ 2 , Im _ 2 ), " (106) 

~ -2 (a + q (i!, lm _ 2 )) ^^.•••.«»-».*»-») (0)a .-(aH- 9 (ix,...,<™- a ))-l/a (1 + Q{z)) 
x— v0+ ' 

^ a .-(a+5(ti,...,* m _a,+ 1 -))-X/2 
k^0+ 



The result Eq (|94| is then still valid at the m th step. 

We still have to verify that this condition is compatible with the regularity of the preceding potentials in the 
chain. For that, note that at the m th step the explicit expression of the resulting potential and its associated 
eigenfunctions are a priori independent of the ordering chosen for the sequence . This means that we can con- 
sider several chains leading to the same extended potential. Consider the set {nx,i\, n m _2,i m — it n m—\i^m— l) 
in which we have regrouped first all the rij associated to a positive ij, that is, we have re-arranged the chain as 

(n'i, i[, —,n' q+ ,i' q+ ,n' q++1 ,i' q++1 , ...,n' m _ 1 ,i' m _-^J = (n^, +,..., n' q+ , +,n q++l , n' m _i, —J, where (n^, n' m _^) is 

a permutation of (m, n m _i) and where q+ = q+ i TO -i). 

Suppose first that q + = m— 1, that is, i'j = +, Vj G {1, m — 1}. Vj G {1, to — 1} we have 1 — q + ^) < 

f Tl % TL % ■ I //-//-/ \ 

and since a > 1, the , ' .', ' r are all nodeless implying the regularity of all the potentials V\ ni ' %1 ''"' n j'*i' in the 
chain. 

Suppose now that q + < to — 1, that is, i' m _i = — ■ The condition a > 1 — q (i 1 ^, i' m _i) = 2 — 
q i' m -2) > 1~<Z ••-> *m-a) implies the absence of node for <^>L 11 l '_ ' m ~ 2 ' m ~ 2 ) and consequently the regularity 

of y«.*i."-,'C-l>«m-l). 

We can iterate this reasoning to go backward along the chain and we see clearly that at each step, the opposite of 
the "total charge" decreases, ensuring the regularity of the associated potential. In this procedure, once we attain the 
level g+, we have eliminated all the negative i'j and retrieve then the first case. 

Consequently, if we satisfy the condition a > 1 — q (ii, i m -i) at the step to, we are ensured of the existence of a 
chain of regular extensions which leads to a regular V^ ni,%1 "'"'' n " m ' lm \ 

IX. EIGENSTATES OF THE EXTENDED POTENTIALS 

We can generalize the calculations made in the two-step case to determine explicit expressions for the eigenstates of 
the successive extensions. For that we adopt the ordering mentioned above and consider the spectrum of the potential 
y(m,+,...,n q+ ,+,n q++1 -,....,n m ,-) where q+ _ q+ (^ ^j^ Using EqflSBJ) and Eq([7HD, we have 
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W< 



(j>m,+ ,-,<f>n q+ ,+,-,4>n g++1 ,-,—,4>n m -,^k I x) (107) 
„ x (m + l) { m/2 +a ) e - { m+l) Z ,2 w ( _ z) ^ ^ ( _ z) ? z -«£-a _^ (z) ? ^ z -a L -« (j8)j £g (z) | x ) , 

that is, with Eq([79]) and Eq(( 



W 







Vii,+> •••7 < r« g+ ,+) •••) 'Png.+l, 



a; (m+l)(m/2+a) e (g + -(m+l)/2) Zz -(m-g + )(Q+m) det ( ^jS >+>•••>««+.+>««+ + 1 , 



where (i, j g {1, m + 1}) 



(ni, + ,..., n q+ , + , n q+ + 1 , — n m ,- 



In the same manner 



1,3 



(108) 



i^-^-z), if j < q+ 
{n 3 -a-i + 2) i _ 1 ^- i+1 L- a - i+1 (z), if g + + 1 < j < m (109) 
i-iy-'L^izUz), if j = m + l. 



that is, 



IE 



('/'t l i,+ 1 -^n !+ ,+ I - I i, + + i-,-^n m - | (HO) 

im (M)/^) eW /J W ( L a ( „ z) ; . f ^ ( _ z) ( z -« e -*L-« +1 (z) , z-" e -^^(z) I x) , 



w 







'«!,+ ) •••) fn q ,+ i <¥n q+ + \,— j Y>n„ 



^ a ,m((m-l)/2+a) z (m- 9+ )(-Q-m+l)) e -(m/2-g+)2 ^ ^(ni,+, 

where e {l,...,m}) 



,n q+ , + ,n q+ + 1 , 



(HI) 



(ni,+, 



,+,n g+ + i, 



(112) 



Eg (1701 gives then 



(«!, + ,. ..,n q+ , + ,n 



?+ + !' 



det * 



(x;w,a) ~ ip (x;oj,a + q(ii,...,i m ))- 



det ^$( ni '+>->"«+'+'™«++i' 



The ^)( ni ' ll ''"'""»* lm ' 1 are then obtained as the product of the modified "gauge factor" (ie fundamental state of 
the initial potential) ipo(x; ui, a + q (ii, i m )), multiplied by the ratio of two polynomials, the polynomial denom- 
inator det 

.,n q+ ,+,n q++ i, ,....,n TO , n being common to all the eigenstates of the considered extension. The 



(113) 



(ni ,ii ,n m 



orthogonality condition on the \p k 
orthogonal family with respect to the weight 



ensure that the det ^> 



[ni,+,...,n g+ ,+,n 



?+ + !> 



constitute an 



r ct+q(ii,...,im) p —z 



det 
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X. SHAPE INVARIANCE OF THE EXTENDED POTENTIALS 

Consider a chain of regular strictly isospectral extensions 

A(v ■ ) W" 1 :* 1 ') Aft/" 1 '? 1 ' n m-K*m-lh 

y 1 ^ yfni,.!) £J a ... nm " m _» V( ni ' <1 »--'*"" <m > ) (114) 

The superpartner of the potential V^ ni ' H ''"' n " m ' t ' n ^(x;ui, a) is given by [22l. l23j 

^(m,i 1 ,...,n m ,i m )j a .. W)0 j = y (ni ' il '-"' n ™' im) (a;;cj,o)+2(4 ni ' <1 ''''' nm ' im) (a;;w,a)) / , (115) 

w^ ll ' ll '"'' nrn ' lrn \x]LJ 1 a) being the RS function associated to the ground level of \Z( n i> l i>---.nm,«m) (£; Q ( w ) = 
Since (see Eq©) 

(n± .ii ,. . . ,n m ,im ) / \ (Wi ,ii .....nm- 1 1 ) / \ Tim A m 

W uc;w,a = —v„ '■ Mx;u;,a) ? — : : 7 ? — : : r , 

^ ni f'-"'^- 1 '^(x; W ,o)-«;^' tl '-"' n -'- 1, **»- l) (x; W ,o) 

(116) 

we have with Eq (|69|) 



_ 2 f gn„,i m (^,a) \ fU8 s 



/ \ ini,ii,...,rt m _i,'i m _i 1 / \ 

(x;cj,a) — Wq (cc;w,a) / 



We proceed by induction. Suppose that the potential V^ ni ' %1 ''"' nm ~ 1 ' %m ~ 1 ' has the same shape invariance properties 
than V 

7 ( " 1 ^--' n " l - 1, ^- l) (x;a;,a) = r (ni ' <1 '-"'^- 1,im - l) (a;;w,a)+2 (tfl^ 1,il,,,,,nm - 1,<m - l) (x;w,o)) / (119) 

_ v( ni ' il '-' nm - 1 ' im - 1 \x;u,a+l) + 2oj. 
As shown in [19j, this is effectively the case for m — 2. It results 

y{mM,..,n m ,i m )^,. fl ) _ F (nx,ii,...,» TO ,i m )^. W; a + i) + 2w - 2 f A< ni,<1 '"" n "» i "'> (x; w, a)) ' , (120) 

where 

A (n ll i 1 ,...,n m)im ) ( ) = g "-^ (CJ ; a) . : 

+lu (n 1 ,i 1 ,...,n m _ 1 ,i m -i) (a; . Wj fl) + ^.ix,...^-!.^-!)^. w? ffl + 1} . 



(121) 



Suppose also that 

A (n 1 ,i 1 ,..., nm _ 1 ,i m _ 1 )( a . ;Wja j = Q) ( 122 ) 

which again is verified for m = 2 fl~oj ] . We can then write 

(ni,il,...,n m _2,im-2) / , i\ (rai,il,---i"m-2,*m-2) / \ /inol 

«„ m -i,i m -i ; (a;;w,a+l) = ; (z;w,a) (123) 

fn m _i,i m _i (^i a) 

(ni,il,...,nTn — 2)im— 2) / \ (ni,21,...,Tl m — 2iim— 2) / \' 
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Moreover we have 

(rii.ii n m _i,i m _i)/ \ (ni.ii,... ,n m -2,im-2) t \ ^m-i-'m-il"'") 



[x;uj,a)—WQ (x;uj,a) 

« (x;u;,a) = — w„ (x:uj,a)-\ — -, — : ■ 1 , 



and (see Eq®) 



i.- ± , - ± , ...,.0 m _ z^m — z / / \ 

(x;uj,a)— . {x;oj,a) 

(124) 



£n m ,i m (u>,a+ 1) - £n m _i,i m _i (w,a+ 1) = £n m ,i m (w, a) - ^7i m _i,t m _i (w,o) . (125) 
Inserting Eq (|124[) . Eq (|123p and Eq (|125|> into Eq (|121|l . we obtain after a little elementary algebra 

A(ni,<i,...,n m ,i ra )( x . WjO ) = 0> ( 12 6) 
By induction, this property is verified for every m. Then Eq (|120j) becomes 

y{n u i^...,n m ,i m ) ^. ^ fl) = j/K,t (3.. W; a + 1) + 2w , (127) 

that is, the potential has also the same shape invariance properties than the isotonic potential. The translational 
shape invariance of the isotonic potential is then hereditary in such chain of extensions. 

XI. CONCLUSION 

We have shown that the "generalized SUSY QM partnership" that we have previously elaborated in a one-step 
scheme can be extended in a multi-step formulation. We have proven the necessary conditions to obtain chains of 
regular extensions of the isotonic potential of arbitrary length and have given explicit expressions for their eigenstates. 
We also established explicitely the hereditary character of the shape invariance properties of the isotonic potential 
which are common to all the potentials in a given chain. 

The case of the generic potentials of the second category of primary TSIP [2l|, namely Poschl- Teller or Scarf 
potentials, can be considered in the same way. This work is in progress and a forthcoming paper is in preparation. 
For the first category exceptional TSIP, due to the strict isospectrality constraint, the only interesting case in which 
we can envisage such chain of extensions is the ERKC potential [2(|. This is the object of further investigations. 
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